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SPHERE BEHAVIOR AND THE MEASUREMENT OF WIND PROFILES 


SUMMARY 


Relatively large, smooth, superpressure, spherical balloons used to 
measure vertical wind profiles in the atmosphere are aerodynamically unstable 
and experience erratic motions when operating in the supercritical range of 
Reynolds numbers. These motions are larger in amplitude than the corresponding 
scales of atmospheric motions. These balloons cannot be used to measure a 
vertical wind profile accurately. The purpose of the present investigation is 
to answer the following questions: (1) What are the characteristics of the aero- 
dynamic forces which cause the spurious balloon motions? (2) Can a smooth, 
superpressure balloon be modified to eliminate the spurious motions and pro- 
vide accurate wind profile measurements? 

Forces associated with aerodynamically unstable balloons were in- 
vestigated by solving the equations of motion and evaluating individual force 
terms using high precision radar data obtained by tracking balloons as they 
ascended through the atmosphere. Aluminized, superpressure balloons made 
of 0. 0127 mm (1/2 mil) Mylar with diameters of 1. 22 m (4 ft) , 2 m (6. 56 ft) , 

2. 13 m (7 ft), and 2.44 m (8 ft) were used. Spring-loaded valves maintained 
a superpressure of 600 to 800 N/m 2 (6 to 8 mbars) , which was sufficient to 
assure sphericity even in the largest expected wind shear conditions. 

As regards smooth, superpressure balloons operating in the super- 
critical range of Reynolds number, the following conclusions were reached: 

1) Statistical properties of the spurious motions are predictable; 2) There were 
no systematic differences observed in the spurious motions associated with 
balloons of different diameters; 3) Lift and drag forces are nearly symmetrical 
along the horizontal axes with a mean near zero; 4) Spurious motions in the x-z 
and y-z planes occur with equal magnitude and frequency, and have a mean in 
each plane of near zero; 5) Distributions of drag coefficients are nearly sym- 
metrical with average values between 0. 25 and 0. 30; 6) The drag force does 
not contain any preferred periods of oscillation; 7) The lift force is predominantly 
horizontal, increases in magnitude with balloon size, and is distributed sym- 
metrically in the x and y directions with means near zero; 8) The lift coefficient 
usually varies between 0. 03 and 0. 05 with some values exceeding 0. 07, while 
mean and RMS values decrease with an increase in Reynolds number in the 
supercritical range; 9) Below a nondimensional wave number of 0. 013 the lift 
force is somewhat organized, but appears random at higher nondimensional 



wave numbers, and; 10) There are no definite periodicities for vortex forma- 
tion and separation, but shorter period changes occur with more regularity than 
longer period changes. 

A smooth, superpressure, spherical balloon was modified to provide 
accurate vertical wind profile measurements. The following conclusions were 
reached: 1) The addition of conical, surface roughness elements 0. 0762 m 
(3 in. ) in base diameter, 0. 0762 m (3 in. ) high and spaced at 0. 19 m (7. 5 in. ) 
apart, and a point mass of 100 g essentially eliminated the spurious motions; 

2) Wind shear does not significantly influence the balloon's response capabilities 
since the maximum expected wind shear is small compared to the ratio of balloon 
ascent rate to balloon diameter; 3) The drag coefficient for the rough balloon is 
almost independent of Reynolds number and has a large value at all altitudes; 
and 4) Wind profile measurements averaged over 25- to 50-m altitude intervals 
can be made with an RMS accuracy of about 0. 5 m/ sec or less to an altitude of 
18 km when the balloon is tracked by the FPS-16 radar. 


INTRODUCTION 

Statement of the Problem 


Balloons are used extensively to measuxe vertical wind velocity profiles 
in the atmosphere. Altitude and wind speed resolution are dependent upon 
balloon properties such as mass and its distribution, shape, size and surface 
roughness, and upon the capability of the system used to track the balloon. 
Relatively large, e.g. , 2-m diameter, smooth, superpressure spheres ex- 
perience spurious oscillations as they rise even through calm air [ 1 , 2] . 
Spurious motions are also observed when these balloons rise through the atmos- 
phere, making it impossible to measure true wind motions on a small scale. 
Measurements of small-scale changes in the horizontal wind along the vertical 
axis are needed in the design and operation of space vehicles, and in meteoro- 
logical problems such as turbulent diffusion, clear air turbulence, energy trans- 
fer, wave motion and forecasting. 

The purpose of this study is to investigate the forces on superpressure, 
spherical balloons which cause spurious oscillations, and to develop an aero- 
dynamically stable, spherical balloon wind sensor capable of providing accurate 
small-scale wind measurements to an altitude of 18 km. 
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Review of Previous Studies 


Reynolds number and drag are commonly used as parameters in in- 
vestigations of sphere motion [3, 4] . The Reynolds number for a sphere is 

defined as ^ ^ — where p is fluid density, V the velocity relative to the sphere, 

M 

d the diameter of the sphere, and p the coefficient of viscosity of the fluid. Drag 
is defined as Jp C D A| V| V where is the drag coefficient and A is the cross- 

sectional area of the sphere. At small Reynolds number, transition and flow 
separation takes place near the equator of a sphere, resulting in a large eddying 
wake. The boundary layer upstream of the point of separation is laminar. As 
Reynolds number is increased, transition from laminar to turbulent flow and 
the separation point move downstream, resulting in reduced wake and drag. 

These facts are fundamental to almost all investigations relating to spheres. 

Bacon and Reid [ 5] summarize previous work on sphere motions and 
drag up to that time. Most investigations had been conducted in wind tunnels, 
and it was found that the type of support such as pendulum, back spindle, wires, 
etc. , influenced greatly the results obtained. Results obtained from hydrogen- 
filled rubber balloons rising in calm air were even more irregular than those 
obtained in tunnels. Drag curves obtained by different investigators were not 
consistent. A satisfactory explanation of the differences had not been given. 

Bacon and Reid investigated the influence of various factors such as tunnel 
turbulence, sphere support, etc. , in an attempt to account for the discrepancies. 
Spheres were dropped in air for the purpose of correlating with tunnel data. 

The spheres were dropped from an altitude of 610 m (2000 ft) during very 
light wind conditions and tracked by theodolites. Among other types, rubber 
balls 0. 305 m (1 ft) in diameter were used. Highly erratic data were obtained 
which were attributed to a patch on the surface of the sphere. The patch caused 
a noticeable increase in drag. Polished wood spheres 30. 5 to 38. 1 cm (12 to 
15 in. ) in diameter were also used; these gave consistent results for Reynolds 
numbers up to 1. 5 x 10 6 . No information was given on the density of these 
spheres. 

Bacon and Reid pointed out that, even though the Reynolds numbers were 
equal for two geometrically similar objects, the Reynolds number of the tunnel 
itself changes with velocity, and it is impossible to get complete dynamic similarity 
because the airflow is not geometrically similar when referred to the dimensions 
of the spheres. It was found that turbulence had a pronounced effect on similarity 
conditions. Consistent results between wind tunnel and falling sphere data were 
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not obtained. In comparing falling- sphere and tunnel data, Bacon and Reid con- 
cluded that tests in free air demonstrated that existing wind tunnels could not 
simulate turbulence conditions prevailing in the atmosphere. 

Hoerner [ 6] further investigated the influence of support methods in wind 
tunnels on the drag of spheres and, in addition, investigated the effects of sur- 
face roughness. He found that surface roughness had a decisive effect on sphere 
turbulence measurements particularly at supercritical Reynolds numbers, and 
that it was necessary to use accurately designed and highly polished spheres in 
making turbulence measurements. He concluded that at supercritical Reynolds 
numbers the degree of surface roughness was far more important than Reynolds 
number effects in determining drag. Hoerner also found that surface roughness 
composed of grains of sand caused the transition (critical) Reynolds number to 
decrease, and that the drag becomes independent of the Reynolds number for 
very rough spheres. 

Platt [7] extended the work on spheres in wind tunnels, using spheres of 
several sizes and wind tunnels of different types. He correlated drag with pres- 
sure variations over the sphere. This work was extended even further by Dryden, 
et al. [8, 9] by consideration of a new variable, viz. , the scale of turbulence. 

Hot wire anemometers and spheres were used to measure turbulence in tunnels. 

It was concluded that the critical Reynolds number of spheres depends on the 
scale of turbulence as well as its intensity. 

Roshko [ 10] studied the wake behind 2-dimensional bluff bodies. He con- 
sidered relationships between the wake and the body including wake scale, fre- 
quencies, energy and interference between wakes. Experiments were conducted 
in a 50. 7- by 50. 7-cm (20 in. ) low-turbulence wind tunnel at the California 
Institute of Technology to investigate coupling between the wake and potential flow 
immediately behind bluff cylinders. He found large pressure fluctuations about 
two plate widths behind plates in the region where vortices form alternately. He 
demonstrated that vortex dynamics had a strong effect on the base pressure be- 
hind a circular cylinder by using a splitter plate to inhibit vortex formation. 
Although Roshko's work contributed greatly to the understanding of the relation- 
ships between bodies and properties of their wakes, empirical relations were 
still required. 

Fung [ 11] studied lift and drag forces on circular cylinders at super- 
critical Reynolds numbers. Experiments were conducted in a tunnel using a 
cylinder approximately 0.305 m (1 ft) in diameter and 1. 83 m (6 ft) long. Lift 
is the force in the direction perpendicular to the velocity vector, and drag is 
parallel to the velocity vector. He observed that the lift and fluctuating part of 
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Lift Coefficient, 


the drag were random phenomena. He found that, in general, the variation from 
the mean did not exceed 7 percent, but on occasions it reached 15 percent. The 
lift coefficient curve for a cylinder taken from Fung's work is shown in Figure 1. 



0 0.2 0.4 0.6 0. 8 1.0 1.2 1.4 

Reynolds Number x 10 ^ 


FIGURE 1. THE AMPLITUDE OF 


Fung performed a power spectrum 
analysis of the lift force which is a ran- 
dom function of time. The power spectra 
varied from sample to sample, but the 
spectra were not found to vary signifi- 
cantly with Reynolds number in the 
supercritical range. There were no 
large peaks observed in the spectra. A 
similar analysis was performed for the 
drag force with peaks appearing in the 

spectra near a Strouhal number ( ^ ; 

n = cycles/ sec, d = diameter of cylinder, 
v = velocity) of 0. 3. Fung concludes by 
stating that the most important feature 
of the forces induced on a circular cy- 
linder by vortex shedding in the super- 
critical Reynolds number range is the 
randomness. He further concludes that 
the mechanism of vortex shedding is 
poorly understood and defies theoretical 
treatment although Roshko's work opens 
up new vistas of research. 


THE LIFT COEFFICIENT FOR A Using expandable neoprene 

FIXED CYLINDER (AFTER FUNG) balloons of different sizes and weights, 

Killen [ 1] performed experiments in a 
large hangar where the wind was calm. Some of the balloons had a weight at- 
tached to simulate a meteorological instrument package usually carried aloft by 
balloons, while others had nothing attached. All the balloons ascended in an 
erratic fashion and without consistency between runs. In some cases, after 
terminal velocity was reached, the balloons appeared to oscillate about the 
vertical axis, while in other cases their distance from the vertical axis increased 
with altitude. There were no systematic motions observed. Killen concluded 
that aerodynamic lift forces were responsible for the observed horizontal motions, 
and that small-scale changes in the horizontal wind field could not be accurately 
measured by balloons. Similar tests by Murrow and Henry [2] in the same or 
similar hangar confirmed Killen' s results. 
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Preukschat [ 12] investigated the motion of rising and falling spheres. 

He found that the wavy motions of a rising sphere became more significant as 
the density ratio (density of sphere to that of the liquid) decreased. The motions 
resulted from a low mass-to-force ratio where the force resulted from non- 
symmetrical vortex shedding. Kising spheres were observed to follow a 2- 
dimensional wavy path which sometimes changed into helical motion. A distinct 
wave length was found for both types of motion and a Strouhal number of 0. 05 
calculated which is about 3 times smaller than that for circular cylinders. For 
low density ratios, the drag coefficient scattered randomly about a mean value, 
but as the density ratio approached 0. 93, the drag coefficient approached the 
steady-state value independently of Reynolds number. 

MacCready and Jex [ 13] performed experiments with spheres rising and 
falling in air and water. They found that sphere motions depend directly on 
Reynolds number, relative mass of the sphere to the displaced fluid (density 
ratio), rotational inertia, surface roughness, sphericity, and random orienta- 
tion. Preukschat's conclusion that the amplitude of the lateral motions in- 
creases as the density ratio decreases was confirmed. Fairly regular zigzag 
or spiral motions were observed at subcritical Reynolds numbers where the 
wake separation is laminar. Wave lengths of about 12 times the diameter of the 
sphere were observed. At supercritical Reynolds numbers where the wake 
separation is turbulent and the wake smaller, an irregular spiral was observed. 
MacCready and Jex concluded that spherical balloons operating in the subcritical 
Reynolds number range, and at higher Reynolds numbers with the addition of 
surface roughness elements, are good wind sensors, although limits of the 
capabilities for such balloons were not given. 

Leviton [ 14] proposed that a 2-m diameter, smooth, superpressure 
sphere be used to make detailed vertical wind profile measurements. The FPS- 
16 high precision radar would be used for tracking the sphere. Leviton con- 
cluded that wind speeds averaged over about 30 m altitude could be measured 
with an RMS accuracy of 0. 3 to 1.0 m/ sec. Scoggins [ 15, 16 ]-, MacCready 
and Jex [13], and Rogers and Camitz [17] investigated the behavior of 
smooth, superpressure balloons and found that it is impossible to make 
accurate wind measurements using such balloons in the supercritical 
Reynolds number range because of their erratic behavior associated with 
vortex shedding. Scoggins [15] proposed a solution by adding large, 
conical, surface roughness elements to control vortex formation and 
separation thus stabilizing the wake, and by attaching a small point mass 
to provide rotational stability and increase the density ratio. Initial experi- 
ments of what is now called the "Jimsphere" were conducted by Scoggins, 
and results later confirmed by Reid [ 18] , MacCready and Jex [ 13] , and 
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Rogers and Camitz [17] . Rogers and Camitz, observing motions of the 
Jimsphere with Doppler radar, found a systematic oscillation with a period 
of 4. 5 sec. This oscillation is almost entirely removed by the data reduc- 
tion procedure employed by Scoggins [ 19] . 

Eckstrom [ 20] developed fabrication techniques and performed a rather 
comprehensive investigation of the Jimsphere and its response capabilities. Drag 
coefficient curves were determined from wind tunnel experiments. The drag 
coefficient was found to be essentially independent of Reynolds number confirming 
calculations by Scoggins [ 16] . Eckstrom determined the theoretical response 
of the Jimsphere and concluded that it was capable of sensing wind motions with 
dimensions of only a few meters. • 


Origin and Importance of the Present Study 


Studies began in 1961 to develop a system for making improved wind 
profile measurements by Jiusto [21] and Figge, et al. [22] . While these 
studies were not entirely successful, another program by Henry, et al. [ 23] , 
employing a smoke producing rocket and cameras, provided wind profile 
measurements of the required accuracy and altitude resolution. The smoke 
trail method turned out to be too expensive (about $3000 per velocity profile) 
for operational use, dependent on good weather conditions, and required ex- 
tensive data reduction procedures. The present study to develop an accurate 
balloon wind sensor began in 1963 after it was realized that the above methods 
were unsatisfactory. Full-scale flight data using balloons of different con- 
figurations were obtained by Scoggins [15] at the NASA Marshall Space Flight 
Center, Huntsville, Alabama, and at the Kennedy Space Center, Florida. Data 
obtained by tracking the balloons with an FPS-16 high precision radar form 
the primary data source for the present investigation. Experimental measure- 
ments by radar were restricted to the Kennedy Space Center. The spherical 
balloon wind sensor developed as a result of the present investigation has been 
in routine use at the two major national test ranges since late 1964. Some of 
the measurements have been published by Scoggins and Susko [ 24] . 

Small-scale motions may cause excitation of dynamic bending modes, 
fuel slosh, and the control system of large space vehicles [ 25] . Geissler [26] 
presented control laws for the Saturn space vehicle which may be used to 
alleviate some of the wind effects; however, atmospheric winds remain per- 
haps the single most important meteorological factor in the design and operation 
of space vehicles. In fact, the establishment of an optimum control law requires 
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a knowledge of the spectral composition of the wind field. Improved wind 
velocity profile measurements are badly needed in various areas of meteoro- 
logical research. Danielsen [ 27] , Weinstein, et al. [ 28] , Endlich [ 29] , and 
others have conducted meteorological research employing small-scale wind 
measurements of the type provided by the wind sensor developed as a result of 
the present study. These measurements are expected to continue to be used in 
fundamental investigations of space vehicle response and meteorological re- 
search. 

EQUATIONS OF MOTION OFA RISING SPHER I CAL BALLOON 

Unstable Balloon 


When a relatively large, smooth, superpressure balloon ascends through 
the lower atmosphere, the Reynolds numbers are in the supercritical range. 

The boundary layer is turbulent, and the separation point is downstream near 
the back of the sphere thus producing a wake with a small diameter. Separation 
does not occur at the same points around the sphere at all times because of 
build-up and decay of vortices or tripping of the boundary layer by surface ir- 
regularities. Presumably, the wake meanders around on the back of the sphere 
producing forces which cause erratic motions as the balloon ascends. In this 
section it is assumed that the balloon moves through a calm atmosphere so that 
all motions of the balloon relative to the air are identical to those measured 
relative to a fixed point on the earth. 

An ascending balloon experiences drag, lift, buoyancy, gravity, and 
pressure forces. The variation of pressure over the surface of an accelerating 
sphere results in a force which acts as an inertia force if the balloon is imparting 
momentum to the air, and acts as an accelerating force when the air is imparting 
momentum to the balloon. The magnitude of this force is given by the product 
of one half the mass of the displaced fluid (apparent mass) and the acceleration 
of the balloon [ 30] . As a sphere moves through an ideal fluid, the momentum 
of the fluid is increased in all directions from the sphere to a distance of 
several sphere diameters, diminishing from the sphere outwards. Essentially 
the same flow field surrounds a sphere moving through a real fluid except that 
a wake forms behind the sphere. Assuming the. wake to be a cylinder extending 
to infinity with a diameter equal to that of the sphere, the momentum computed 
from ideal fluid theory in this cylinder was found to be approximately 4 to 5 per- 
cent of the total momentum imparted to the fluid. This is the largest wake 
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expected; therefore, the apparent mass may be treated as constant. Moreover, 
because of symmetry, the apparent mass associated with a sphere is a scalar 
quantity. In a calm atmosphere the balloon imparts momentum to the air so 
that the apparent mass term acts to increase the inertia of the balloon and may 
be considered as part of the mass of the system. 

From Newton's second law of motion 

(M + M a ) V B = SF = B+ L+D, (1) 

After substitution for the forces, this equation becomes 

(M+M a )V b = g[ VOL(p g -p a )+M B ]+ ip a C L A|V| (Vxn)+ip a C D A| V|V, 

(2) 

where M = mass of the balloon and gas 
M a = apparent mass 

= mass of balloon 

VOL = volume of balloon 
A = cross-sectional area of balloon 
P a = density of air' 

p = density of gas 

s 

V = velocity of air relative to the balloon 

V = acceleration of the balloon 

= drag coefficient 
C = lift coefficient 

Li 

g = acceleration of gravity 

n = a unit vector perpendicular to V and in the plane of the lift force. 
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The first term on the right-hand side of equation (2) contains buoyancy and 
gravity; however, these forces will be treated jointly. These forces always act 
in the vertical direction. The second and third terms on the right-hand side of 
equation (2), lift and drag forces , respectively, act in perpendicular directions 
with drag being defined as having a direction along the relative velocity vector. 
The lift force may assume any direction in the plane perpendicular to the drag 
vector. 


Equation (2) may be written in component form for a three-dimensional 
orthogonal coordinate system. The coordinate system used in this report as- 
sumes a flat earth with x positive toward the east, y positive toward the north, 
and z positive upward. In this coordinate system the component equations of 
equation (2) are: 

x-component 

(M+M a )V Bx = £p a C L A|V|(Vxn) *7+ £p & C d A|V|V x (3) 

y-component 

(M+ M A )V By = iP a C L A|V| (Vx7) * T+ sP a C D A IVI V, (4) 

z -component 

ip a C L A m (Vxn)*k 

(5) 

j, and k are unit vectors. 

The component forces are considered positive when they point in the 
positive direction along the axes. The vector lift force is perpendicular to the 
vector drag which restricts it to a particular plane at a given instant of time, 
but does not restrict its direction within that plane. Thus, the direction of the 
component lift forces and, therefore, that of the vector lift force as well, must 
be such that equations (2) through (5) balance. 

The equations of motion for an unstable balloon rising through a changing 
wind field will not be considered here. Balloon motions caused by unstable lift 
forces are generally an order of magnitude larger than true small-scale atmos- 
pheric motions (over altitude layers of less than approximately 300 m) and, 
therefore, it is not possible to measure these motions with an unstable balloon. 


(M+M a )V Bz = g[ VOL(p-p a )+M B ] 


+ K C D A|VIV z 


Along the x, y, aqd z axes, respectively, i, 
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Stable Balloon 


A stable balloon is by definition one that rises vertically along a straight 
line in a calm atmosphere. The lift forces are zero, and the only forces acting 
on the balloon are buoyancy (defined to include gravity) and drag. The vector 
equation of motion is 

(M + M a )V b = l[VOL(p g -p a )+ Mg] + iP a C d A|V|V (6) 

which may also be written as 

(M+M A )V Bz =g[VOL(p g -p a )+M B ]+ip a C D A|V|V z , (7) 

since there is no motion in the horizontal directions. These equations are valid 
only for a stable balloon rising through a calm atmosphere in which case the 
balloon imparts momentum to the atmosphere, and the apparent mass term acts 
to increase the inertia of the system. 

We now consider a stable balloon rising through a changing wind field. 

The vector equation of motion is 

mv b = 7[VOL (p g - p a )+ m b ]+ ip a c D A| V|V - m a V. (8) 

The apparent mass term has been separated from the inertia term since mo- 
mentum may be transferred from the balloon to the air or vice versa. If mo- 
mentum is transferred from the balloon to the air, the last term in equation (8) 
will act as an inertia term, but if momentum is transferred from the air to the 
balloon, this term acts to accelerate the balloon. 

OBSERVED AERODYNAMIC CHARACTERISTICS OF SMOOTH SUPER- 
PRESSURE BALLOON WIND SENSORS RISING FREELY IN THE 

ATMOS PHERE 

Description of the Wind Sensors 

In meteorology, as well as in the design of space vehicles, it is desirable 
to obtain accurate wind measurements to an altitude of at least 15 to 20 km and, 
in some cases, to a much higher altitude. In order to obtain wind measurements 
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to a given (maximum) altitude, the balloon wind sensor required is a function 
of its buoyancy-to-weight ratio. For a sphere made of a given material, the 
larger the diameter, the higher will be its floating altitude. However, it is more 
difficult to launch larger balloons especially during high winds. It is generally 
better to use a balloon no larger than required to reach the desired altitude. 

A 2 -m diameter, superpressure sphere made of 5 mil, aluminized mylar 
with two release valves and one fill valve made of plastic, weighs approximately 
300 grams and will float at an altitude in excess of 20 km when inflated with 
helium. A sphere of this type is shown in Figure 2. The pressure release 

valves are spring-loaded and maintain 
a superpressure of approximately 600 to 
800 N/m 2 (6 to 8 mbars). The balloon 
is used as a passive target and there is 
no instrumentation attached. 

Observed Spurious Motions and 
Their Dependence on Reynolds 
Number 


A series of measurements made 
at approximately 1 -hour intervals with 
a 2 -m diameter sphere of the type de- 
scribed above as the wind sensor is 
shown in Figure 3. The balloons were 
tracked with the FPS-16 radar. These 
measurements were made at Cape 
Kennedy, Florida, on January 3, 1961. 
The data reduction method used to com- 
pute wind velocity from radar position 
data is given in Appendix A. The solid 
line superimposed on the dotted profile 
represents the wind profile measured by the rawinsonde system. This profile 
represents winds averaged over approximately 600 m and is shown for com- 
parison only. 

At high Reynolds numbers the boundary layer associated with the balloon 
is turbulent while at low Reynolds numbers the boundary layer is laminar. 
Reynolds numbers are called supercritical in the turbulent flow regime, and 
subcritical in the laminar flow regime. At supercritical Reynolds numbers, 
flow separation takes place on a small area near the back of the sphere resulting 
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Altitude (km) 



Wind Speed (m/sec) 


FIGURE 3. SCALAR WIND SPEED PROFILES MEASURED AT 
APPROXIMATELY 1-HOUR INTERVALS USING SMOOTH, 2-m DIAMETER 
SUPERPRESSURE BALLOON (SOLID LINE IS RAWINSONDE MEASURED 

PROFILE) 


in a small wake. As vortices build up and separate, the separation region 
moves around on the back of the sphere, causing the wake to tilt. Tilting of the 
wake, in turn, causes lift forces which are responsible for the observed 
spurious motions. At subcritical Reynolds numbers, flow separation takes 
place farther upstream (near the sphered equator) , resulting in a large wake. 

In this case the formation and separation of vortices do not cause the wake to 
tilt significantly, so the lift forces are much smaller than the drag forces. Thus, 
at subcritical Reynolds numbers a sphere is aerodynamic ally more stable than 
at supercritical Reynolds numbers. 
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In Figure 3 the motions of the balloon are observed to be erratic to an 
altitude of approximately 12 km, and above this altitude the erratic behavior 
disappears. It is at this altitude that the flow regime for the 2-m diameter 
sphere used experiences transition from a supercritical (turbulent boundary 
layer) to a subcritical (laminar boundary layer) flow regime. Coincident with 
this transition, the turbulent wake associated with the sphere increases in size 
and becomes more stable. Thus, at subcritical Reynolds numbers the lift forces 
are drastically reduced in magnitude, the drag force becomes larger and more 
stable, and the balloon becomes a better wind sensor. 

EVALUATION OF THE EQUATIONS OF MOTION FOR 
SMOOTH SUPER PRES SURE BALLOONS FROM EXPERIMENTAL DATA 

Solution of the Equations 

A solution of equations (2) through (5) and evaluation from experimental 
data provide an improved understanding of the forces which act upon a rising, 
aerodynamic ally unstable spherical balloon and which produce spurious or self- 
induced motions. The spurious motions were determined by defining the "true" 
motions of the balloon by use of a numerical filter function (see next section) , 
then subtracting these motions from the measured motions. Since the spurious 
balloon motions are several times larger than atmospheric motions on the same 
scale as shown from measurements by independent methods, the definition of the 
filter function is not too critical. The filter function could be replaced by a 
moving arithmetic average without degrading the results. 

With the assumption of no vertical air motion, the motion of the balloon 
relative to the air is given by the vector sum of the ascent rate of the balloon 
and the horizontal components of the spurious motions. Spurious motions in the 
vertical direction are small compared to the ascent rate and are neglected. Ac- 
celerations of the balloon were calculated from the measured speeds for the x, 
y and z directions. Accelerations computed from the spurious balloon motions 
defined above were almost identical to those computed directly from the meas- 
ured speeds. This is because acceleration caused by spurious balloon motions 
are an order of magnitude larger than those caused by a change of true wind 
speeds which the balloon experiences as it rises. Using the measured accel- 
erations, the sum of the forces acting along each axis was determined and then 
resolved into individual force components. 
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Density, p, apparent mass, M^, buoyancy, B, and the sum of the forces 
along each axis, F^, F , and F^, were computed as follows: 
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where r is radius of the sphere, T is temperature, is molecular weight, R* 

is the universal gas constant, and p is pressure. Pressure and temperature are 
obtained from radiosonde flights. 

The coordinate system used and 
an illustration of the resolution of the 
velocity vector into its components are 
shown in Figure 4. Except for signs, 
the direction of the measured velocity 
vector may be replaced by the direction 
of the drag force vector since the drag 
vector is defined as being opposite the 
measured velocity vector. The fol- 
lowing relationships are obtained di- 
rectly from Figure 4: 

V H V 

cosc^p^j-, cos e =^-. (15) 

FIGURE 4. COORDINATE 
AXIS SHOWING RESOLUTION OF 
WIND VECTOR AND THE RELATIONSHIP 
(16) BETWEEN FORCE AND WIND VECTORS 


Solving for V^ we get 

V x = | V | cos (p cos 6 , 
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where 

i 

| V| = [ V 2 + V 2 + V 2 ] 2 . 
x y z 

The expression V may be replaced by -D and V by D . Since cp and 6 are de- 

X X 

termined by the component wind speeds, equation (15) may be substituted into 
equation (16) giving 


D = — 1 D | 
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cos 4> cos 6= — | D | ( 
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Similarly, equations for the drag force components in the y and z directions are 
given by 

_ V 

D = - |D| — (18) 

y i vi 

and 

V 

D = - |D| ~=± . (19) 

z I V| 


The magnitude of the drag force, | D| , is given by 

— ► — i 

| D 1 = | F| cos a = [ D 2 + D 2 + D 2 ] 2 

1 x y z J 


( 20 ) 


where 


cos a 


F V+F V+F V 
xx y y z z 

|F| | V| 


( 21 ) 


and 

i 

| F | = [ F 2 + F 2 + F 2 ] 2 . 
x y z 

The term a is defined as the angle between F and V, and is a function of the 
magnitude of the lift force (see Figure 4) . 
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The lift force components are given by 


L = 

F - D 

(22) 
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F - D 

(23) 
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F - D 
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z z z 

and the magnitude of the lift force, | L| , and the horizontal component of the lift 
force, L r , by 

|L| = [L 2 + L 2 + L 2 ]® (25) 

1 x y z 

and 

L h = [L* + L 2 y f . (26) 

The direction of the horizontal lift force, \I> , is given by 

l L I 

\ 1 ! = tan -1 |^j + quadrant correction, (27) 

I xl 


where the quadrant correction is determined by the signs of Ly and L^. 

The drag and lift coefficients, C and C , respectively, are given by 

D Li 


C D * 


2 1 D I 


7T p r 2 I V| 2 
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and 


C L = 


2|L| 

7T p a r 2 | V| 2 


(29) 
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Characteristics of the Forces 


Analysis of many ascents of superpressure, spherical balloons which 
have a smooth surface and are aerodynamically unstable at supercritical 
Reynolds numbers, show that only the statistics (means, variances, etc. ) of 
self-induced or spurious motions are predictable. Spiral, planar, and random 
motions have been observed from balloons rising in the atmosphere as well as 
from spheres rising or falling in water, but with no regularity [ 13] . The types 
of motion observed appear to be a function of the formation and separation of 
vortices which are random events. Vortex formation and shedding phenomena 
associated with flow around cylinders show similar results in the supercritical 
flow regime [ 11] . Because of the random nature of the spurious motions and 
thus of the forces associated with aerodynamically unstable balloons, statistical 
methods (empirical probability distributions and spectra) are used to investigate 
the behavior and interrelationships of the forces acting on the rising balloons. 

The average ascent rate of balloons rising in the atmosphere in the 

supercritical flow regime is essentially constant. The balloons are tracked in 

space as a function of time, and velocities computed as a function of altitude 

(Appendix A) . Because of the relatively constant ascent rate, it is possible to 

relate space and time through the relationship V At = AZ. For example, wave 

z 

number spectra may be converted to frequency spectra and vice versa. 

In the analysis which follows, balloons with diameters of 1. 22 m (4 ft), 

2 m (6. 56 ft) , 2. 14 m (7 ft) , and 2. 44 m (8 ft) are used. These are designated 
as tests 1761, 1762, 1767, and 1768, respectively. All measurements were 
made on the same day at Cape Kennedy, Florida. The same FPS-16 radar was 
used to track the balloons. The balloons were all made of 0. 0127 mm (J mil) , 
aluminized Mylar and inflated with helium. The superpressure of 600 to 800 
N/m 2 (6 to 8 mbars) is sufficient to maintain sphericity even in large shear 
layers. The 1. 22-m (4-ft) diameter sphere goes from supercritical to sub- 
critical Reynolds numbers at an altitude near 5 km, and the larger balloons near 
13 km. Spurious motions are analyzed below these altitudes for the respective 
balloons. 

Because of the symmetry of a spherical balloon, statistical properties 
of the erratic motions in the horizontal plane should be independent of direction. 
This was verified along the x (x-z plane) and y(y-z plane) axes. In the 
analysis that follows, data are given only for the y-axis. Wind shears are an 
order of magnitude smaller than V /d and would not be expected to significantly 
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influence the spurious motions of the balloons. Superpressure balloons of 
different sizes rising through the atmosphere give similar results in the super- 
critical flow regime. This is because of Reynolds number similarity, which 
means that characteristics of the flow over the balloons are independent of 
viscosity. 

Composite Forces. Composite forces represent the sum of all forces 
acting on the balloon along a given axis. These forces, which include the effects 
of apparent mass, were computed using equations (12) through (14). The 
magnitude of the inertia term in equation (14) is small compared with buoyancy; 
therefore, the composite force in the z-direction is essentially equal in mag- 
nitude to buoyancy, which is treated in a later section. 

Probability distributions of the composite forces along the y-axis are 
shown in Figure 5 for the four balloons. All forces have been normalized by 


Percent 



FIGURE 5. PROBABILITY DISTRIBUTIONS OF NORMALIZED 
COMPOSITE FORCES ALONG y-AXIS 
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| F| . Cumulative frequency distributions or ogives corresponding to these dis- 
tributions and plotted on probability paper are shown in Figure 6. The distribu- 
tions are normal or Gaussian as indicated by the straight lines and the 95 percent 



FIGURE 6. CUMULATIVE FREQUENCY DISTRIBUTIONS (OR OGIVES) 

OF NORMALIZED COMPOSITE FORCES ALONG y-AXIS 

confidence limits. The method for computing confidence limits for a cumulative 
frequency distribution, assuming the theoretical distribution is normal (Gaussian) , 
is given by Hald [ 31] . The 95 percent confidence limits were computed at ap- 
proximately the mean, ± 1 standard deviation, and ± 2 standard deviations for 
each cumulative frequency distribution. These limits were then transferred to 
the probability distributions via the scale of the variable since it is common to 
both figures. This is not a common statistical procedure, but it does provide 
an indication of how much the probability distributions can vary and still remain 
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within the 95 percent confidence limits. The 95 percent confidence limits are 
shown on all probability distributions. 

The probability distributions of the composite forces are nearly sym- 
metrical with a mean near zero. This indicates that the spurious motions, which 
primarily determine the accelerations over small altitude intervals, are also 
symmetrical; i. e. , the balloon oscillates back and forth about its mean path as 
it rises through the atmosphere. 

Individual F o rce Terms' 

Buoyancy . Buoyancy is defined as the net upward force due to the 
difference in densities between the balloon system (including gas) and the air. 

It is convenient to consider gravity in this term since it acts in the vertical and 
does not appear in any of the other terms. At lower altitudes, buoyancy is a 
large force which decreases with altitude as the air density decreases. Floating 
altitude is reached when the buoyancy force becomes zero. 

Dra g. Probability distributions of the y-components of the drag 
force given by equation (18) are shown in Figure 7. Probability distributions 


Percent 



FIGURE 7. PROBABILITY DISTRIBUTIONS OF NORMALIZED 
DRAG FORCES ALONG y-AXIS 


of the forces in the x-direction were similar to those in the y-direction. The 
vertical component of the drag force decreases in magnitude with altitude and 
becomes zero at the floating altitude. The magnitude of the vertical component 
of drag is an order of magnitude greater than the magnitudes of the horizontal 
components. This is caused by the greater component of velocity relative to 
the balloon in the z-direction. 

Figure 8 shows probability distributions of the drag coefficients. The 
distributions are reasonably symmetrical with the average value of the drag 



coefficient varying between approximately 0.32 for the 2-m (6. 56-ft) diameter 
sphere to 0. 24 for the 2. 13- and 2. 44-m (7- and 8-ft) diameter spheres. This 
variation in average drag is caused by variation in balloon size and ascent rate. 
Average values of the drag coefficient for a 2-m (6. 56-ft) diameter sphere cal- 
culated by Scoggins [ 16] agree very closely with the most probable value of 0. 32 
obtained above for supercritical flow. The most probable values for the other 
spheres are about a factor of 2 larger than those determined in wind tunnels [ 3] . 

Spectra of the drag force and the drag coefficient are shown in Figures 9 
and 10. Spectra of the drag force do not show any preferred periods of oscillation 
except at low Strouhal numbers, a fact which is caused by a large steady-state 
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Normalized spectral density 


value. In the case of an ascending 
balloon, the Strouhal number is given by 
the ratio of balloon diameter to the 

d 

vertical wave length, — . This is a non- 

A 

dimensional wave number which is in- 
dependent of the ascent rate of the balloon. 
For all tests except 1767, nondimensional 
wave numbers above 0. 01 contain about 
equal energy. This means that oscilla- 
tions in drag over small ascent distances 
occur randomly. Spectra of the drag co- 
efficients for the larger balloons (tests 
1767 and 1768) show peaks near non- 
dimensional wave numbers of 0. 013 and 
0. 030. The overall trend in all spectra 
is a decrease in energy with an increase 
in nondimensional wave number. 

According to Blackman and Tukey 
[ 32] , the accuracy or statistical con- 
Nondimensionai wave number, i x io 2 fidence of a computed spectral estimate 

depends upon the length of record and 
the number of lags used. Spectral 

FIGURE 9. SPECTRA OF estimates are considered as a chi-square 

DRAG FORCES variate with the number of degrees of 

freedom, k, given approximately by 

2T /T , where T is the length of record and T is the maximum lag. In the 
n m n m 

spectrum analysis which follows, the number of degrees of freedom is approximately 

15. According to Blackman and Tukey, the computed spectral estimates for 15 

degrees of freedom will be within approximately ± 50 percent of their value with 

95 percent confidence. These confidence limits are shown on all spectra. 

Lift. Probability distributions of lift in the y-direction are shown in 

Figure 11. Distributions in the x-direction were similar to those in the y-direction. 

Magnitudes of the lift forces in all directions increase with balloon size with the 

magnitude in the z-direction being several times smaller than in the x- and y- 

directions. This is caused by the orientation of the wake. The z-component of 

lift is always negative because of the large force of buoyancy which causes a 

large composite force, F , in the negative z-direction, which, in turn, determines 

z 
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FIGURE 10. SPECTRA OF 
DRAG COEFFICIENTS 


the sign of F . Probability distributions 

of lift in the z-direction are shown in 
Figure 12. For all balloons, the 
magnitudes of the component lift forces 
in the x- and y-directions are almost sym- 
metrical about zero and have the same 
distribution in both directions, indicating 
no statistical preference for direction. 

Probability distributions of the 
magnitude of the lift vector are shown in 
Figure 13. Distributions of the horizontal 
component of lift were very similar to 
those of the total lift force indicating that 
the lift force is predominantly horizontal. 
The average magnitude of the lift force 
increases with balloon diameter. 

Distributions of the lift coefficient 
defined by equation (29) are shown in 
Figure 14. The most frequently occurring 
value varies between approximately 0. 03 


Percent 




FIGURE 11. PROBABILITY DISTRIBUTIONS OF NORMALIZED LIFT ALONG y-AXIS 
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Percent 


and 0. 05; however, the distributions for 
tests 1761 and 1768 are flat with a 
relatively high percentage of occurrences 
of coefficients as large as 0. 07. 

Curves of the mean and RMS 
values of the lift coefficient versus 
Reynolds number in the supercritical 
range are shown in Figure 15. The 
mean and RMS values of the lift coef- 
ficients , which differ very little in mag- 
nitude, decrease with an increase in 
Reynolds number. Fung [ 11] obtained a 
similar relationship for circular cy- 
linders at Reynolds numbers below ap- 
proximately 4 x 10 5 (see Figure 1). The 
curves in Figure 15 have been extended 
to lower Reynolds numbers to indicate 
the sharp decrease in the lift coefficient 
which occurs at subcritical Reynolds 
numbers. 


Spectra of the total lift force are 
shown in Figure 16. Spectra of the 
horizontal lift force were almost identical 
to those of the total lift force. Spectra of 
the total lift force show relatively constant 
energy at nondimensional wave numbers 
above 0. 01, but do show a tendency to- 
ward increased energy at lower non- 
dimensional wave numbers. This indicates 
that changes in the lift force at lower non- 
dimensional wave numbers are somewhat 
more organized than at higher nondimen- 
sional wave numbers. Stated differently, 
the lift force experiences somewhat organized large-scale variations with super- 
imposed random variations. Spectra of the horizontal lift force show the same 
features indicating the lift force is primarily horizontal. 



|F| 


FIGURE 12. PROBABILITY 
DISTRIBUTIONS OF NORMALIZED 
LIFT ALONG z-AXIS 


Spectra of the lift coefficient are given in Figure 17. Low energy is 
shown for tests 1762 and 1767 at nondimensional wave numbers of 0. 013 and 
0. 035, respectively, with an overall decrease in energy with an increase in 
nondimensional wave number. 
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FIGURE 15. MEAN AND RMS VALUES OF THE LIFT COEFFICIENT 
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FIGURE 16. SPECTRA OF TOTAL LIFT FORCES 

Spectra of the changes in direction of the horizontal lift force defined by 
equation (27) are shown in Figure 18. With the exception of test 1767, which has 
a peak near a nondimensional wave number of 0. 030, the only significant feature 
of these spectra is an overall increase in energy with an increase in nondimensional 
wave number. This means that short period changes in the direction of the 
horizontal lift force occur with more regularity than the longer period changes. 


Relationship Between Total and Drag Forces 


It was assumed previously that the drag vector acted in the direction of 
the relative velocity vector, but there was no restriction placed on the direction 
of the total force vector. The angle between these two vectors is a function of 
the magnitude and direction of the lift force. Only a small percentage of the 
angles is greater than 16 degrees (cosine a = 0. 97) , with the percentage in- 
creasing sharply for the smaller angles. These large angles are indicative 
of an unstable balloon and are caused by large lift forces. 
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FIGURE 17. SPECTRA OF LIFT COEFFICIENTS 
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Nondimensional wave number, — x 10 

FIGURE 18. SPECTRA OF CHANGES IN HORIZONTAL 
DIRECTION OF LIFT FORCES 


THE EXPERIMENTAL DESIGN OF AN ACCURATE 
SPHERICAL BALLOON WIND SENSOR 


Based on work done by Killen [ 1] , Scoggins [ 15, 16] , MacCready and 
Jex [ 13] and Rogers [ 17] , it became apparent that the instability of a sphere 
moving through a fluid is associated with an unstable wake causing erratic lift 
and drag forces. Also, it is evident from the results presented in the preceding 
section that a smooth, superpressure balloon large enough to reach a 20-km 
altitude is not a satisfactory wind sensor for the troposphere (surface to 12 km) . 
Elimination of the lift forces is required in order to stabilize the wake. Although 
much has been learned about the aerodynamic behavior of such balloons, a com- 
plete understanding still does not exist, and it is not possible to describe ana- 
lytically the spurious motions observed when the balloon rises through the atmos- 
phere. Even if an exact analytical description were available, the interaction 
between the wind field and the balloon would have to be known before accurate 
wind measurements could be made. This would require a prior knowledge of the 
wind field which is the very thing we desire to measure. 

With this background an approach was taken, based on present knowledge 
and sound theoretical principles, to conduct an experimental program aimed at 
developing an accurate spherical balloon wind sensor. The problem appears to 
be that of controlling vortex formation and separation which, in turn, controls 
lift and drag. To do this, balloons with varying types of roughness elements 
were fabricated and flown, and the data were analyzed to determine what in- 
fluence the roughness elements had on the spurious motions. Roughness ele- 
ments change the velocity profile in the boundary layer which causes earlier sep- 
aration as in the case of a subcritical sphere. Based on results obtained from 
each series of tests, additional configurations were flown and data analyzed. The 
process was repeated until the desired results were obtained. Eight configura- 
tions were tested. 

Wind profiles were measured by the smoke trail method [ 23] and used as 
a reference in evaluating the effectiveness of roughness elements to eliminate 
spurious motions. An experimental analysis (Appendix B) shows that the smoke 
trail wind data has an RMS error of generally less than 0. 5 m/ sec. 

In addition to smoke trail wind profile data, the 1. 22-m (4-ft) diameter, 
smooth, superpressure balloon was also used as a reference in each series. 

Even though this balloon experiences self-induced motions [17, 33, 34] , it can 
be used as an indicator of performance when operating in the subcritical flow 
regime for comparison with other balloons. The 1. 22-m (4-ft) diameter sphere 
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goes from supercritical to subcritical Reynolds numbers at an altitude of ap- 
proximately 4 to 6 km. Below this altitude the balloon operates in the super- 
critical Reynolds number regime and is unstable. An appreciation for the dif- 
ference in a balloon operating at subcritical Reynolds numbers and at supercritical 
Reynolds numbers may be gained by comparing the spectra for the 1. 22-m 
( 4— ft) diameter balloon below the transition altitude with that above. For ex- 
ample, spectra of the small-scale scalar balloon motions (see subsection 
following, "Analysis of Measured Balloon Motions," for definition) are shown 
in Figure 19 between 2. 5 and 6 km (supercritical) and between 6 and 13. 3 km 
(subcritical) . 

In the development of an improved wind sensor, emphasis must be placed 
on the characteristics of the measured small-scale motions. In this section the 



FIGURE 19. SPECTRA OF 
SMALL-SCALE SCALAR MOTIONS 
FOR SAME BALLOON IN SUB- AND 
SUPERCRITICAL FLOW REGIMES 


jffects of various sizes and shapes of 
roughness elements on the variability and 
spectra of small-scale motions are pre- 
sented. The effects of roughness elements 
and attitude orientation on aerodynamic 
stability are discussed. Drag curves for 
various balloon configurations and a com- 
parison of smoke trail and balloon- 
measured wind profile data are also 
discussed. 

The Experimental Program 

The experimental program con- 
sisted of flying spheres of different con- 
figurations in sequence on the same day. 

The types of balloons were alternated to 
provide comparisons between as many 
balloon configurations as possible. The 
test series were conducted when the radar 
was available rather than during chosen 
weather conditions. It turned out that 
low wind speeds existed during all the 
test series. All balloons were launched 
from the same location and tracked by 
the same radar. Pre- and post-calibrations 
of the radar were accomplished. 
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Data for three series of balloon flights are shown in Table I. A number 
of additional series similar to these were also flown. Table I is self-explanatory 
except for the column headed "type. " This column refers to the configuration; 
e. g. , 2m, 3T(or F) , 9S refers to a balloon 2 m in diameter (2m) with 3-inch 
(0. 076-m) high truncated or full cones (3T or F) spaced 9 inches (0. 2286 m) 
apart on the sphere (9S) . The "weight" column refers to the balloon alone and 
does not include the gas. 


TABLE I, DATA FOR SERIAL BALLOON FLIGHTS 


Test Number 

Release Time (Z) 

Type- 

Weight (g) 


Series 1 - 

June 9, 1964 


4378 

1800 

4 ft smooth 

163 

4379 

1916 

2m, 3T, 9S 

301 

4380 

2045 

2m, 3F, 7S 

301 

4381 

2210 

4 ft smooth 

164 

4382 

2345 

2m, 4F, 7S 

310 

4383 

0115 

2m, 3T, 9S 

308 

4387 

1900 

Smoke trail 



Series 2 - J 

une 11, 1964 


4384 

1400 

4 ft smooth 

163 

4385 

1521 

2m, 3F, 7. 5S 

302 

4386 

1640 

2m, 4F, 8S 

308 

4393 

1806 

4 ft smooth 

165 


Series 3 - June 

12 and 13, 1964 


4405 

1700 

4 ft smooth 

164 

4497 

1830 

2m, 2T, 6S 

400 

4498 

1950 

2m, 2T, 6S 

300 

4499 

2245 

4 ft smooth 

164 

4500 

0016 

2m, 2T, 6S 

600 

4502 

0130 

2m, 2T, 6S 

1000 

4503 

1900 

Smoke trail 



* These units have not been converted to S. I. units. 
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Analysis of Measured Balloon Motions 


Definition of Small-Scale Motions . Since the problem is to develop a wind 
sensor that will respond accurately to small-scale changes in the wind field, 
emphasis must be placed on analyzing the higher frequency components of the 
measured profiles. For this section small-scale motions are defined as those 
which are not included in the rawinsonde measured profiles. This definition is 
not of too much importance as long as the longer wavelengths included in the 
definition of small-scale motions are longer than the longest wavelength of the 
self-induced motions. With this definition, wavelengths as long as approxi- 
mately 500 m are included, extending well beyond the wavelengths of all spurious 
balloon motions. 

The small-scale components of motion were separated from the profiles 
by defining a filter function by the relationship 

S R (f) = |F| 2 S g (f) (30) 

where S is the spectrum of the wind profile measured by the smoke trail method 

and assumed to be the true wind speeds, S is the spectrum of the wind profile 

measured by the rawinsonde system, and F is the filter function. The filter 
function, and the weights which must be applied at 50-m intervals to a detailed 
(25-m interval) profile in order to produce a smooth profile containing the same 
frequency content as would have been measured by the rawinsonde system, is 
shown in Figure 20. The weights were obtained by taking the cosine transform 
of the filter function, then normalizing. The small-scale motions (wavelengths 
less than approximately 500 m) are obtained by applying the filter weights in a 
step-wise manner to a detailed wind profile to define a smooth profile, then sub- 
tracting this smooth profile from the detailed profile. The higher frequency 
motions represented by the difference between the smoothed and detailed pro- 
files are used to assess the stability of the various balloons. Equation (30) , 
with smoke trail data replacing rawinsonde data and balloon data replacing 
smoke trail data, was used to define erratic balloon motions used in the pre- 
ceding section. 

The resulting filter function was similar to that shown in Figure 20 but 
with a cutoff near a wavelength of 300 m. The same weights are used but applied 
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FIGURE 20. FILTER FUNCTION 
USED FOR DEFINING SMALL-SCALE 
MOTIONS 


at 25-m intervals rather than 50-m 
intervals. 

General Observed Features from 
Serial Ascents . Profiles of balloon 
motions for the test series presented in 
Table I will not be shown here. Most 
of the profiles exhibited features some- 
what similar to those shown in the 
preceding section but with much less 
spurious motions. The profiles meas- 
ured with the rough balloon configurations 
show variable scatter over limited altitude 
bands, although the scatter was in most 
cases not sufficient for detection when 
plotted. For this reason it was necessary 
to use quantitative methods in assessing 
performance of the balloons. Remarkable 
consistency was noted between adjacent 
profiles indicating that true larger scales 
of the wind motions were being measured. 


Performance of the V arious Balloon Configuration s. The performance of 
each balloon was evaluated by computing the variance and spectra of the small- 
scale motions for each flight and comparing with data obtained with the 1. 22-m 
(4-ft) diameter sphere and smoke trail. The energy or variance associated with 
each test shown in Table I and the altitudes over which each spectrum was com- 
puted, are shown in Table II. Variations in the amount of energy contained in 
the small-scale motions from one profile to another are clearly shown in Table II. 


The variances associated with each profile using the rough balloon con- 
figurations were generally larger than those for the 1. 22-m (4-ft) diameter 
smooth balloon in the subcritical Reynolds number flow regime and for the smoke 
trail profiles. From visual observations of the ascent of several of the balloons, 
it was noted that they were experiencing some rotation. It was further observed 
that when rotations occurred the balloon appeared to experience somewhat larger 
lateral motions. The roughness elements had eliminated a large percentage of 
the spurious motions by forcing flow separation and reducing lift, but the rotation 
apparently produced a small lift force which still had to be eliminated. This was 
done by attaching a small mass at a point on the sphere which displaced the center 
of gravity downward, thus providing a stabilizing torque. This also increased the 
density ratio which has a stabilizing effect [ 12] . Spectra for Jimsphere type 
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TABLE II. VARIANCE AND ALTITUDE DATA CORRESPONDING TO 
BALLOON FLIGHTS PRESENTED IN TABLE I 


Test 

Number 

Altitude 
Range (km) 

Zonal (V ) 
' x ' 

Variances 

Meridional (V ) 
y 

Scalar (|V|) 







S 

eries 1 - June 9, 

1964 


4378 

3. 0-14. 0 

0. 246* 

0. 211* 

0. 211* 

4379 

2. 0-20. 0 

0.432 

0. 435 

0.381 

4380 

3. 5-20. 0 

0.431 

0.442 

0. 340 

4381 

2. 0-8. 0 

0. 634 

0.715 

0.732 


8.0-15. 0 

0. 159 

0. 246 

0. 106 

4382 

8. 0-15. 0 

0.208 

0. 289 

0. 210 


3. 5-20. 0 

0.494 

0. 601 

0. 541 

4383 

2. 0-20. 0 

0.489 

0. 587 

0.457 


& 

5ries 2 - June 11, 

1964 


4384 

2. 0-6. 0 

0. 576 

0.462 

0. 522 


6. 0-13. 3 

0. 095 

0. 186 

0. 145 

4385 

2. 0-16. 0 

0. 514 

0.419 

0.594 

4386 

3.0-18. 0 

0.450 

0. 452 

0. 452 

4393 

1. 8-5.5 

0.420 

0. 682 

0. 334 


6. 5-14. 0 

0. 173 

0. 178 

0. 144 


Serie 

;s 3 - June 12 anc 

1 13, 1964 


4405 

5.0-10. 0 

0. 517 

0. 482 

0. 502 

4497 

6. 0-14. 0 

0. 162 

0. 204 

0. 226 

4498 

6. 0-14. 0 

0. 204 

0. 243 

0. 226 

4499 

2,5-7. 0 

0. 366 

0. 531 

0. 378 


7. 0- 14. 0 

0.217 

0. 254 

0. 134 

4500 

2. 5 - 15. 0 

0. 185 

0. 193 

0. 184 

4502 

3. 5 - 13. 0 

0. 161 

0. 160 

0. 158 


* m 2 /sec 2 
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2m, 3F, 7. 5S with a 100-g mass attached, smoke trail, and a 1. 22-m (4-ft) 
diameter sphere operating at subcritical Reynolds numbers were almost iden- 
tical. This means that RMS errors caused by erratic motions associated with 
Jimsphere balloon type 2m, 3F, 7. 5S with a mass of 100 g attached are no larger 
than those associated with the smoke trail and 1. 22-m (4-ft) diameter sphere 
data. As noted by Rogers and Camitz [ 17] this Jim sphere configuration experi- 
ences a regular spiral motion with a wavelength near 25 m. This spiral motion 
is eliminated by the data reduction method presented in Appendix A [ 17] . 

Drag Coefficient Curves 


The value of the drag coefficient is of primary interest in computing the 
response of a balloon to a change in wind. The magnitude of the drag force, 
which accelerates the balloon in the presence of a wind shear, is directly pro- 
portional to the magnitude of the drag coefficient. Thus large drag coefficients 
are desirable. 

The magnitude of the drag coefficient increases with the size of the wake. 
The drag coefficient measured in wind tunnels for smooth spheres is small at 
supercritical Reynolds numbers where the wake is small, increases sharply in 
the transition range, and has a much larger and almost constant value at sub- 
critical Reynolds numbers where the wake is large [3,4]. A minimum is ob- 
served at supercritical Reynolds numbers immediately above the transition 
range. The addition of roughness elements to a spherical balloon causes separa- 
tion near the equator, resulting in a large wake and a large drag coefficient at 
all. Reynolds numbers. 

Drag coefficient curves are shown in Figure 21 for Jimsphere type 2m, 

3F, 7. 5S calculated from full-scale flight data and measured in a wind tunnel 
using a scaled model. For comparison, a curve measured in a wind tunnel using 
a smooth sphere is also shown. The drag coefficient calculated from full-scale 
data is larger at all Reynolds numbers than was determined in the wind tunnel. 

Preukschat [ 12] found that the ratio of the density of a sphere to the 
density of the fluid was an important factor in determining how the sphere will 
move through the fluid. A sphere constrained in a wind tunnel has an effective 
infinite density, while a sphere such as type 2m, 3F, 7. 5S rising in the atmos- 
phere has an average density ratio of about 0. 4. This difference in density ratios 
may account, at least in part, for the difference in the drag coefficient curves. 
Another factor which may account for some of the difference is that a free-rising 
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Drag coefficient, 



FIGURE 21. DRAG COEFFICIENT 
CURVES FOR JIMSPHERE TYPE 
2m, 3F, 7. 5S AND FOR A SMOOTH 
SPHERE IN A WIND TUNNEL 



FIGURE 22. JIMSPHERE TYPE 
2m, 3F, 7. 5S 


balloon can adjust very quickly to changes 
in forces caused by a nonsymmetrical or 
tilted wake, thus keeping the drag vector 
pointing opposite the direction of motion 
of the balloon. Scoggins [ 16] investigated 
the drag coefficient for a rough balloon 
with a density ratio of about 0. 8, and 
computed a drag curve almost identical 
to the curve determined in a wind tunnel 
with a model of type 2m, 3F, 7. 5S 
(Fig. 21). This supports the hypothesis 
that the difference in the curves in 
Figure 21 for the rough spheres is ac- 
counted for, at least in part, by the 
difference in density ratios. 

Dynamic Response Characteristics 
of Aerodynamically Stable Rough 
Balloons to Vertical Wind Shears 

Eckstrom [ 20] , using the drag 
coefficient curve determined in wind 
tunnel tests and presented in Figure 21, 
investigated theoretically the response of 
Jimsphere type 2m, 3F, 7. 5S to wind 
shears. This balloon is shown in Figure 
22. Eckstrom found that the Jimsphere 
would respond with a lag of less than 
1 m/ sec to the maximum expected wind 
shear of 0. 2 sec -1 at all altitudes below 
13 km. Below 10 km the lag is less than 
0.3 m/sec. Using the drag curve de- 
termined from full-scale flight data, these 
lag errors would be even less. In either 
case, the balloon responds rapidly enough 
to permit the measurement of wind shears 
over a few balloon diameters. In practice, 
when radar tracking errors and data re- 
duction are considered, average winds 
over 25 to 50 m can be measured with an 
RMS error of generally less than 0. 5 
m/ sec. 
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A typical scalar wind speed profile measured with the Jimsphere is 
shown in Figure 23. Spurious balloon motions associated with smooth, 
superpressure balloons of the same size as the Jimsphere, and shown 
in Figure 3, have been eliminated. The same data reduction method 
was used in both cases. 



FIGURE 23. SCALAR WIND SPEED PROFILE MEASURED 

WITH JIMSPHERE 



CONCLUSIONS 


The following specific results were obtained in this study: 

1. The statistical properties of spurious motions associated with smooth, 
superpressure, spherical balloons operating at supercritical Reynolds numbers 
are predictable. 

2. There were no systematic differences observed in the spurious 
motions associated with balloons of 1. 22 m (4 ft) , 2 m (6. 56 ft) , 2. 13 m (7 ft) , 
and 2. 44 m (8 ft) in diameter. 

3. Composite, lift, and drag forces along the horizontal axes are nearly 
symmetrical with a mean near zero. 

4. Spurious balloon motions in the x-z and y-z planes occur with equal 
magnitude and frequency, and have a mean in each plane of near zero. 

5. Distributions of the drag coefficients for balloons of various diameters 
are nearly symmetrical with the average value varying between 0. 30 for the 

1. 22-m (4-ft) diameter sphere to 0. 25 for 2. 13-m (7-ft) and 2.44-m (8-ft) 
diameter spheres. 

6. The drag force does not contain any preferred periods of oscillations. 

7. The lift force is predominantly horizontal. 

8. Magnitudes of the component lift forces in the x and y directions are 
almost symmetrical with a mean near zero. 

9. Magnitudes of the component lift forces increase with balloon size 
with the magnitude in the z-direction being several times smaller than in the x- 
and y-directions. 

10. The most frequently occurring values of the lift coefficient vary 
between 0. 03 and 0. 05 with the value reaching as high as 0. 07 ten percent or 
more of the time. 

11. Mean and RMS values of the lift coefficient decrease with an in- 
crease in Reynolds number in the supercritical range. 
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12. The lift force experiences somewhat organized large-scale varia- 
tions with superimposed random variations. While there are no preferred modes 
of variation, there is a tendency toward organized variations at nondimens ional 
wave numbers below 0. 013. 

13. Assuming that changes in direction of the horizontal component of 
lift are indicative of vortex formation and separation, there are no preferred 
periods for these events; however, the shorter period changes occur with more 
regularity than the longer period changes. 

14. The angle between the total force and drag vectors usually does not 
exceed 15 degrees. 

15. The addition of conical roughness elements 0. 076 m (3 in. ) in 
height, 0. 076 m (3 in. ) in base diameter, and spaced 0. 19 m (7. 5 in. ) apart, 
and a point mass of 100 g to a spherical balloon produces stabilization and 
essentially eliminated spurious motions. The roughness elements produce 
more and smaller vortices than those associated with smooth balloons. The 
superposition of the smaller vortices onto the larger ones produces an equal 
distribution of energy in the lift forces over all frequencies. 

16. Damping and inertia of the rough balloon (Jimsphere) with a point 
mass of 100 g attached prevent significant rotation and thus eliminate lift at higher 
nondimensional wave numbers. 

17. The maximum expected wind shear of 0.2 sec -1 is small compared 

with V /d, where V is the ascent rate of the Jimsphere and d its diameter; 
z z 

therefore, wind shear would not significantly affect the balloon’s response 
capabilities. 

18. The drag coefficient for the rough balloon is almost independent of 
Reynolds number and has a large value at all altitudes. 

19. Jimsphere type 2m, 3 F, 7. 5S with a mass of 100 g attached provides 
wind profile measurements averaged over 25- to 50-m altitude intervals with an 
RMS accuracy of about 0. 5 m/sec or less to an altitude of 18 km when tracked by 
the FPS-16 radar. 
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Suggestions for Further Research 


Although this investigation has provided much information on sphere 
behavior and the forces producing this behavior in the supercritical Reynolds 
number flow regime, there are still several fundamental problems to be solved. 
Of particular interest are the phenomena of the formation and separation of 
individual vortices, more explicit relationships between drag, lift, buoyancy, 
and apparent mass effects, the influence of extreme wind shears and turbulence 
on the response of a sphere, and the influence of the wake on the motions of 
the sphere. Further research is also desired in developing a fast-rising wind 
sensor that can reach 18 km within a few minutes (or a system which will other- 
wise measure a vertical wind profile in a short time) rather than the 1 hour re- 
quired by the present configuration. Finally, data handling procedures need 
refining to improve the quality of the measured wind data. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, March 6, 1967 
160-44-04-00-62 
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APPENDIX A 


DATA REDUCTION PROCEDURE FOR SUPERPRESSURE SPHERICAL 
BALLOONS WITH DISCUSSION ON ACCURACY OF WIND DATA 

Edit and Computation Procedure 

The processing of radar data to get wind data is a very important 
function for any system used. Accuracy of the wind data often depends upon 
the methods employed in processing the basic tracking data which, in turn, 
are often limited by available computing facilities. The procedure for pro- 
cessing fast rate (0. 1-sec interval) tracking data presented here was developed 
for use on an IBM 7094 computer [ 19] . While in all probability it is not an 
optimum program in terms of providing the most accurate wind data possible, 
it does give satisfactory results for most applications and requires only a 
nominal amount of computer time. There are essentially two parts to the 
procedure; one is editing the radar tracking data to eliminate erroneous points, 
and the other is processing the edited data to get winds. Each of these aspects 
of the program will be covered separately. 

Data Edit Procedure . The data edit procedure is quite simple. The 
problem is to eliminate values which are not accounted for by trend or random 
error, and to replace them by values which are reasonable. The first step is to 
fit the measured values of r, 0 , and 4 > (slant range, elevation and azimuth re- 
spectively) to a 9-point, first degree equation by the method of least squares. 

If the residual at the mid-point of the 9-point interval is greater than 3 times 
the assumed RMS error (provided by the radar manufacturer or obtained by 
experiment) , the point is replaced by the average value over the interval. In 
order to prevent including stray points in the computation of the mean, the next 
point immediately following the 9 points used in the curve fit is compared with 
the 9-point mean, and if the difference is 5 times greater than the assumed 
RMS error, the point is replaced by the mean. This procedure is repeated for 
each 0. 1-sec step until all the data have been edited. 

Wind Computation Procedure . The procedure for computing winds using 
the edited radar tracking data is given below as a series of steps. 

a. Convert slant range from yards to meters, and azimuth and eleva- 
tion angles from degrees to radians. 
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b. Compute x, y, and z for each point: 


7 r 

x = r cos 0 cos (<j> - ~ ) (positive east) 

a 

z = r sin 0 (positive vertical) 

y = r cos 0 cos 0 (positive north) 

c. Correct position coordinates x, y, and z for earth's curvature: 


, -1 x 

x = r tan 1 ; — 

c o r + z 

o 


z = [x 2 + (z + r ) 2 + y 2 ! 2 - r 
c 1 ' o' J 1 o 


y = r tan 
J c o 


-l y_. 


[ x 2 + (z + r ) 2 ] 2 


r^ = 6, 373,334 m = radius of earth at Eastern Test Range, Florida. 


d. Determine altitude and corresponding time for 25-m intervals as 
follows: 

1) Find the first z^ which is greater than the desired altitude. 

2) Use this z as the mid-point of a 41-point (4 sec) least squares 
first degree curve fit. 

3) Determine desired y from polynomial. 

4) Repeat 1) , 2) , and 3) for x and y . 

c c 

5) Interpolate linearly between times corresponding to altitudes 
below and above z to obtain time corresponding to z. 

e. Compute component wind speeds x , and y , the scalar wind speed 

n n 

v^, and vertical velocity of the balloon z^, using the following equations: 
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. 

X - X 

n n - 1 

(m/ sec) 

x — 
n 

At 

(n, n - 1) 


<<! 

0 

l 

P 

t 

(m/ sec) 

y n = 

At (n, n- 1) 


Z n “ Z n - 1 

(m/sec) 

z ~ 

n 

At (n, n-1) 

v = 
n 

t"n + 

(m/ sec) 


The computed wind speeds are associated with the top of the layer designated 
by n. 


f. Compute direction from which the wind is blowing, as follows: 
x 

^ = tan" 1 7 — + quadrant correction 

n y 

°n 

The quadrant correction is determined from the sign of x and y as follows: 

n n 


x + 
n 


360 - * 


n 


n 


x + 
n 


n 


+ 180 
n 


x 

n 

. j 180 - * 

y n + ‘ 


n 


y n - 


* 


n 
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Accuracy of Computed Wind Data 


Error analyses have been conducted for determining the accuracy of wind 
data obtained by tracking a balloon by the FPS-16 radar. The accuracy is a 
function of the condition of the radar, the operator’s experience, and the data 
reduction procedure used. A practical approach has been followed here to cir- 
cumvent difficulties in trying to evaluate the errors contributed from each of 
these factors. The approach was to track the same balloon with 2 radars simul- 
taneously and independently, and then to compare the computed wind profiles. 
Three pairs of profiles, one from the Western Test Range (WTR) , and two from 
the Eastern Test Range (ETR) were obtained. A 2-m diameter smooth balloon 
was used at WTR while the Jimsphere was used at ETR. The smooth balloon 
operates in the supercritical Reynolds number flow regime to an altitude of ap- 
proximately 13 km; above this altitude the balloon operates in the subcritical 
Reynolds number flow regime. At subcritical Reynolds numbers, the balloon is 
aerodynamically more stable and, presumably, the true atmospheric motions are 
measured. Thus up to approximately 13 km, the balloon oscillated back and 
forth in both planes (probably a combination of zigzag and irregular spiral 
motions) , causing the radar to "search. " 

The RMS errors in the computed balloon motions were determined by 
assuming that each radar had the same RMS accuracy. Independence is also a 
necessary assumption, but this is obviously satisfied since the radars operate 
completely independently. The RMS error, cr 12 , of the sum or difference of two 
variables (in this case balloon motion or wind speed) having RMS errors of 
a i and <r 2 is given by the relation [ 35] 

ct i ,2 = o \ + °2 • 

Assuming that both radars have the same RMS tracking accuracy, the above 
equation becomes 

cr = 0.707 cr« , 

R 1 >‘ s 

where a is the RMS error in balloon motions caused by the tracking accuracy 

X\ 

of either radar. 

The last equation was used to calculate RMS accuracies of the wind data 
over 2-km intervals for the three profile measurements. The results are shown 
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in Tables A-I, A -II, and A -HI. The smooth balloon (WTR) was aerodynamically 
unstable to an altitude of approximately 12 km, while the Jimsphere (ETR) was 
aerodynamically stable at all altitudes. Thus, below approximately 12 km the 
"search" of the radars is more pronounced in the WTR case than in the ETR 
cases. 


As shown in Table A-I the accuracy of the radars used at WTR improved 
with altitude and is a factor of 2 to 3 times better between 13 and 20 km than 
between 2 and 4 km. In contrast, as shown in Tables A-II and A-III, there is 
only a factor of 1 to 2 improvement in the accuracy of the radars used at ETR 
for the same altitudes when tracking the aerodynamically stable Jimsphere. This 
improvement in accuracy is apparently caused by a decrease in spurious motions 
which causes less searching. At altitudes above 12 km, where both balloons are 
aerodynamically stable, no significant differences were observed. 


TABLE A-I. AVERAGE VALUES OF RADAR COORDINATES AND RMS 
ERRORS IN WIND DATA AS FUNCTION OF ALTITUDE FOR WTR TEST 3002, 

1730Z , MARCH 16, 1963 


ALTITUDE 

INTERVALS 


RMS ERRORS 


AVERAGE VALUES 

V 

x 

V 

y 

(m/sec) 

V 

V 

z 

0 & r 

(m) 

(m/ sec) 

(m/sec) 

(m/sec) 

(deg) (deg) (m) 

2000 - 3975 

0. 49 

0. 66 

0.57 

0. 18 

89.4 32.4 6100 

4000 - 5975 

0.40 

0.71 

0.49 

0. 15 

79.6 25.0 12,900 

6000 - 7975 

0. 40 

0.35 

0.40 

0. 14 

75.2 20.4 21,800 

8000 - 9975 

0. 29 

0.28 

0.27 

0. 17 

73.5 18.4 31,000 

10,000 - 11,975 

'i' 


-I' 

* 

-e 

'r -i' 

12,000 - 13,975 

0. 23 

0.24 

0.24 

0. 23 

74.2 15.3 53,200 

14, 000 - 15, 975 

0. 24 

0.18 

0.24 

0. 29 

73.6 14.1 65,300 

16,000 - 17,975 

0. 16 

0.23 

0. 17 

0. 26 

73.2 12.9 80,900 


* No Data 
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TABLE A -II. AVERAGE VALUES OF RADAR COORDINATES AND RMS 
ERRORS IN WIND DATA AS FUNCTION OF ALTITUDE FOR ETR TEST 
8920, 1402Z, DECEMBER 23, 1964 


ALTITUDE 

INTERVALS 

(m) 

RMS ERRORS 

AVERAGE VALUES 

V 

x 

V 

V 

V 

V 

Z 

4> 

6 

r 

(m/ sec) 

J 

(m/sec) 

(m/sec) 

(m/ sec) 

(deg) 

(deg) 

(m) 

2000 - 3975 

0. 30 

0.39 

0.39 

0. 32 

113. 0 

30. 0 

6450 

4000 - 5975 

0. 23 

0. 16 

0.25 

0. 12 

108. 0 

27. 3 

11,750 

6000 - 7975 

0, 19 

0. 19 

0.23 

0. 12 

114. 7 

25. 0 

17,950 

8000 - 9975 

0. 18 

0.28 

0. 30 

0. 11 

121. 9 

23. 7 

24,475 

10, 000 - 11,975 

0, 19 

0. 18 

0.22 

0. 12 

130. 6 

21. 7 

32,400 

12, 000 - 13, 975 

0. 27 

0. 18 

0. 32 

0. 13 

129. 0 

18. 9 

43, 500 

14, 000 - 15, 975 

0. 32 

0. 18 

0. 32 

0. 15 

125. 1 

17. 0 

55, 300 

16, 000 - 17, 975 

0. 21 

1 

0. 17 

0. 24 

0. 23 

123.4 

16. 3 

65, 100 


TABLE A-HI. AVERAGE VALUES OF RADAR COORDINATES* AND RMS 
ERRORS IN WIND DATA AS FUNCTION OF ALTITUDE FOR ETR TEST 
8940, 1600Z, DECEMBER 23, 1964 


ALTITUDE 

INTERVALS 

(m) 


RMS ERRORS 


V 

X 

(m/ sec) 

V 

y 

(m/ sec) 

V 

(m/ sec) 

V 

z 

(m/ sec) 

2000 - 3975 

0. 30 

0. 23 

0.26 

0. 17 

4000 - 5975 

0. 16 

0. 10 

0.21 

0. 14 

6000 - 7975 

0. 17 

0. 19 

0. 22 

0. 13 

8000 - 9975 

0. 23 

0.27 

0.33 

0. 20 

10, 000 - 11, 975 

0. 16 

0. 19 

0. 21 

0. 13 

12,000 - 13,975 

0. 25 

0. 17 

0. 27 

0. 15 

14, 000 - 15, 975 

0. 28 

0. 14 

0.27 

0. 19 

16,000 - 17,975 

0. 20 

0. 18 

1 

0.20 

0. 18 


* As shown in Table A -II 
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APPENDIX B 

ACCURACY OF SMOKE TRAIL WIND DATA 


Wind data are obtained by the smoke trail method by establishing a 
vertical column of smoke using a small rocket, then determining movement of 
the smoke trail photographically in time. The smoke trail is assumed to move 
with the speed of the wind. 

Analytical procedures for computing wind speeds (horizontal displace- 
ment of the smoke trail during some convenient time period) are presented by 
Houston [ 36] and Foster [ 37] . These procedures are not difficult mathematically, 
but it is difficult to establish analytically the accuracy of the computed wind 
data because of the subjectivity involved. For example, locating the center of 
the trail is a matter of judgement. Errors in wind data are also caused by film 
shrinkage, survey errors, camera pointing inaccuracies, multiple solutions, 
and round-off errors. None of these errors are known precisely. Thus, an 
experimental approach seems to be appropriate to determine the accuracy of 
the wind data. 

The RMS errors in wind data measured by the smoke trail method were 
investigated by Camp and Scoggins [ 38] . They considered two wind profile meas- 
urements from the same smoke trail but used different camera pairs. Table 
B-I was extracted from their work. Errors associated with the profiles were 
determined from the differences between two profiles measured from the same 
trail using different camera pairs. Errors of individual wind speed measure- 
ments sometimes exceed 1 m/sec, but the RMS errors are believed to be 
generally below 0. 5 m/sec as shown in Table B-I. 

Errors in smoke trail wind data caused by vertical air motions have not 
been determined. Identification of points on the trail, which would be required 
for the measurement of vertical motions, has not been possible. It appears, 
however, that features of the smoke trail caused by large wind shears, thin 
layers of relatively fast wind speeds, etc. , generally remain at about the same 
elevation, indicating that vertical air motions in these regions are small. 

Vertical air motions are believed to influence significantly only a very small 
percentage of the wind data measured by the smoke trail method. 
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TABLE B-I. RMS ERRORS (m/sec) IN SMOKE TRAIL WIND DATA 
OVER ALTITUDE INTERVAL 5-13 km FOR THREE PROFILES 


Profile Number 

V 

X 

V 

y 

V 

i 

0. 17 

0. 32 

0. 27 

2 

0, 27 

0. 61 

0.60 

3 

0. 16 

0.43 

0.27 
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